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In this work we show that the standard method to obtain nucleation rate-predictions with the aid of atomistic 
Monte-Carlo simulations leads to nucleation rate predictions that deviate 3 — 5 orders of magnitude from the 
recent brute-force molecular dynamics simulations [J. Diemand, R. Angelil, K. K. Tanaka, and H. Tanaka, 
J. Chem. Phys. 139, 074309 (2013)] conducted in the experimental accessible supersaturation regime for 
Lennard-Jones argon. We argue that this is due to the truncated state space literature mostly relies on, 
where the number of atoms in a nucleus is considered the only relevant order parameter. We here formulate 
the nonequilibrium statistical mechanics of nucleation in an extended state space, where the internal energy 
and momentum of the nuclei is additionally incorporated. We show that the extended model explains the 
lack in agreement between the molecular dynamics simulations by Diemand et al. and the truncated state 
space. We demonstrate additional benefits of using the extended state space; in particular, the definition of 
a nucleus temperature arrises very naturally and can be shown without further approximation to obey the 
fluctuation law of McGraw and Laviolette. In addition, we illustrate that our theory conveniently allows to 
extend existing theories to richer sets of order parameters. 
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I. INTRODUCTION 

Nucleation is the initial step of a phase trans¬ 
formation, a process ubiquitous from the smallest 
to the largest physical length-scales: the fabri¬ 
cation of new nanoscale material for various fu¬ 
ture applications such as the design of new mi¬ 
crochips requires a control on the morphology of 
nucleated materiali:^. On the micrometer scale, 
nucleation is for instance relevant in the devel¬ 
opment of new drugs; here, different nucleation 
conditions can lead to different packing arrange¬ 
ments of the molecular substance, which would 
be pharmaceutically classified as different drugs. 
Therefore a precise handle on the nucleation pro¬ 
cess is deeply required^. On the higher end of 
length scales, research has been conducted- to 
study cloud seeding - the injection of a substance 
into clouds to induce rainfall - offering the poten¬ 
tial to control the weather and the atmosphere. 
On cosmic length-scales nucleation leads to cre¬ 
ation of interstellar dust after super novae. The 
dust modifies the spectrum of interstellar radia¬ 
tion, and its composition and time-scale of forma¬ 
tion are of key interest^. In all these phenomena, 
small nuclei have to overcome an activation bar¬ 
rier to initiate the phase-transformation and to 
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grow deterministically. The rate at which small 
nuclei are capable to escape from the initial fluc¬ 
tuation dominated growth regime is the funda¬ 
mental quantity of interest in nucleation theory. 
This nucleation rate has been predicted through 
a variety of methods ranging from brute force 
molecular dynamics (MD) simulations^to phe¬ 
nomenological, semi-phenomelological and statis¬ 
tical mechanics approaches'^— among others. 

MD simulates all involved atoms in the nucle¬ 
ation process by solving their equations of mo¬ 
tion and can be used to obtain the exact nu¬ 
cleation rates under various initial conditions of 
the metastable phase. Since MD conducts com¬ 
puter experiments with systems of nano sized 
dimensions, up to limits of a few micrometers, 
these studies favour regimes in which the nucle¬ 
ation rate is large, so that despite of the mi¬ 
croscopic system size critically sized nuclei pop 
up frequently enough and adequate statistics on 
the obtained data can be ensured. In the pro¬ 
totypical studied nucleation of liquid drops in a 
metastable gas, MD works in the large super¬ 
saturation regime where the metastable phase is 
very dense. This is in strong contrast to real 
experiments. Most laboratory experiments^Sr— 
are performed at low supersaturations and mea¬ 
sure nucleation rates less than 10^°cm“^s“^. For 
the nucleation of argon the latest experimen¬ 
tal setups through the supersonic Nozzle (SSN) 
experimenf— achieved to study nucleation rates 
up to 10^^cm“^s“^. Only recently, in 2013, large- 
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scale MD simulations^ were capable to operate in 
the experimentally accessible regime set by the 
SSN. In Fig. [T] we have reported the studied nu¬ 
cleation rate in MD simulations of argon in some 
of the literature throughout the past 20 years 
- unless sophisticated large-scale simulations are 
exploited, MD would typically study nucleation 
regimes that differ by 5—10 orders of magnitude 
from the SSN experiment and 12 — 17 orders of 
magnitude from typical experiments. MD cannot 



FIG. 1. Illustration of accessible nucleation rate Jmd in se¬ 
lected molecular dynamics studies^^ of argon from 1997 — 
2013. We marked the first letter of the family name of the 
first author of the study together with the number of used 
atoms. The shaded area shows the experimental accessible 
regime. 

be used to study realistic experimental regimes 
unless tremendeous computational effort is made, 
which is the main drawback of the technique. 
Coarse-grained models allow to bridge the time- 
and length scale gap between atomistic dynamics 
and nucleation as perceived on a coarser scale. 
Such models offer the possibility to obtain pre¬ 
dictions in the experimental accessible nucleation 
regime and offer hence an appealing alternative 
to brute-force MD. 

The development of coarse-grained models to 
capture the relevant macroscale physics of nucle¬ 
ation received tremendeous amount of interest. 
These studies mainly focus on obtaining predic¬ 
tions for the shape and hight of the nucleation 
landscape Vs(nc) that prohibits small nuclei to 
reach the critical size n*. The barrier height 
AVs of Vs is particularily important and enters 
the nucleation rate through an exponential fac¬ 
tor Relatively minor errors in AVs are 

therefore known to give rise to lack in predic¬ 
tions of the nucleation rate of orders of magni¬ 
tude. In the literature k-sTVsiric) is known as the 
Gibbs free energy of formation of a nucleus of 


size Tic in a metastable phase at temperatur T and 
pressure P. A long history of phenomenological, 
semi-phenomenological and statistical mechanics 
models have been obtained to predict the work 
of formation. In particular, the classical nucle¬ 
ation theory^, the internally consistent classi¬ 
cal nucleation theory^, the semiphenomenolog- 
ical model (SP)^^— , the revised SP modeli^, the 
mean-field kinetic theory (MFKT) ^^i^^ , the ex¬ 
tended modified liquid drop dynamic nucleation 
thoerysi^ and many more. In our work we aim to 
obtain quantitative predictions of Vs from first 
principle atomistic Monte-Carlo (MC) simula¬ 
tions thereby avoiding the often crude approxi¬ 
mations that mesoscopic theories exploit to ob¬ 
tain the nucleation landscape Vs- 

Standard MC simulations^^— in the literature 
exploit small atomistic systems with N atoms, 
pressure P and temperature T to sample the 
probability (j>s{nc) of occurence of a nucleus of 
size ric- The choice of the isobaric-isothermal en¬ 
semble in computer simulations mimics the con¬ 
ditions of many real experiments which control 
exactly these variables to induce the nucleation 
process. The link of the microscopically ob¬ 
tained probability 4)s{nc) with the macro-physics 
manifests itself in the relationship Vs{nc) = 
— \og{(j)s{nc))- The probability (j)s{nc) is vanish¬ 
ingly small around the critical size n*, but the MC 
simulations can be biased through importance 
sampling techniques such as umbrella-sampling to 
access regions of the atomistic phase space that 
contain a nucleus of size ric- MC simulations ob¬ 
tain (j)sinc) even in the experimentally relevant 
nucleation regime thereby offering a tool that is 
potentially superior to brute-force MD. In order 
to obtain nucleation rate predictions J, (j)s must 
however enter an adequate evolution equation ob¬ 
tained through supplementary dynamic consider¬ 
ations of the macroscale nucleation process. The 
Zeldovich equationi^ 


df{nc) _ d ^ \ df{nc) 

dt duc [ dric 


fin-c) 


dVsjnc) 

dric 


( 1 ) 


recognized as thermodynamically admissible^ is 
probably the most popular to capture the dynam¬ 
ics of the distribution function /(ric) of a nucleus 
of interest. Here, /(ric) denotes the probability 
that the considered nucleus has size ric and 
is the diffusion coefficient of the diffusion process 
described by Eq. {[T]). The Zeldovich equation can 
be solved approximately in the stationary limit to 
obtain the nucleation rate as J = where 

the prefactor Jg is related to kinetic coefficients 
and the barrier shapei^. 

We will illustrate that when the nucleation 
rate of argon droplets in a supersaturated ar¬ 
gon gas is obtained through the standard formula 
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J = , where AVg is constructed from the 

aforementioned MC technique, a nucleation rate 
is predicted that deviates 3 — 5 orders of mag¬ 
nitude from exact brute-force MD predictions by 
Diemand and coworkers^. This is particularity se¬ 
vere since the simulations conducted by Diemand 
et al. consider the practically relevant snpersat- 
uration regime. Furthermore, the argon parti¬ 
cles are described by the Lennard-Jones model 
and this is one of the most basic and often used 
microscopic models to study the nucleation pro¬ 
cess. The standard approach to obtain the nu¬ 
cleation barrier by MC simulations is therefore 
incomplete. 

We argne that this is an artefact of relying on 
the fully truncated state space = (ric) to de¬ 
scribe the nnclei degree of freedom on the coarse¬ 
grained scale. Onr attempt to improve nucleation 
rate predictions obtained from MC simulations 
has its fnndament in 1966 in the classical work^i by 
Feder et al. who already noticed the need to ex¬ 
tend the aforementioned truncated state space in 
order to account for certain dynamic effects that 
cannot be captured by the use of a single state 
variable. In particular, Feder et al. realized that 
in the nucleation process a significant amount of 
latent heat is released indicating that the nnclei 
exhibit rather strong energy fluctuations on top 
of their size fluctuations. They showed that a 
proper inclnsion of the internal energy Uc into 
the nnclei state space leads to corrections to 
the nucleation rate of 1 — 3 orders of magnitude. 
The nucleation landscape in their work is there¬ 
fore extended to two dimensions such that Vs = 
Vs{nc,Uc)- This is confirmed by other work^^i^^ — 
focnsing on nonisothermal corrections. Another 
argument in favour of equipping the state space 
with additional degrees of freedom traces back to 
the “translational-rotational paradox” initiated 
by Lothe and Pound aronnd 19622S. They real¬ 
ized that nnclei are not static objects, bnt can 
move and rotate within the metastable phase. In 
their semi-quantnm mechanical treatment they 
tried to incorporate the translational and rota¬ 
tional contributions into Vs{nc) withont extend¬ 
ing the nuclei state space and ended up with cor¬ 
rections of 10^^ to the nucleation rate. After in¬ 
corporating their correction, theoretical predic¬ 
tions of the nncleation rate deviated from exper¬ 
imentally obtained rates so outreageonsly that 
the term “translational-rotational paradox” was 
coined. Many attempts^Srii tried to resolve the 
paradox. In this respect the work^ of Reguera 
and Rubi presents an elegant solution by incor¬ 
porating additional degrees of freedom that char¬ 
acterize the nuclei. In particular, they argued 
that a proper inclusion of the momentnm Pc and 
angular momentum Ic in the state space of the 


nncleus leads to much smaller corrections of the 
order of 10^ — 10® for typical nucleation regimes. 
The landscape would hence be extended to Vs = 

Vs(^c; Pcj ^c)' 

The idea of extending the state space to cap¬ 
ture relevant effects is not new and for cer¬ 
tain order parameters, the extended nucleation 
landscape Vs has even been studied by statis¬ 
tical mechanics^—. But in this paper we aim 
to lay the basis for the nonequilibrinm statisti¬ 
cal mechanics of nncleation with extended state 
spaces. We here restrict the study to the spe¬ 
cific choice Xc = («c,Pc,Wc) and to homogeneous 
nucleation in a one-component system. We de¬ 
velop a coarse-graining procedure that is reli¬ 
ably predicting nucleation rates in the low nu¬ 
cleation rate regime where critically sized nuclei 
are sufficiently large so that their dynamics can 
be captnred in the continuum limit through the 
use of Fokker-Planck equations. This Fokker- 
Planck equation involves diffusion in a;c-space 
and describes the dynamics of f{xc), the prob¬ 
ability that at a certain time a specific nucleus 
has property Xc- Our coarse-graining proce¬ 
dure relies on the GENERIC (General Equa¬ 
tion for Non-Equilibrium Reversible-Irreversible 
Coupling) projection-operator technique^i^ and 
yields atomistic expressions for the building 
blocks of the aforementioned Eokker-Planck equa¬ 
tion. In particnlar, we obtain a fundamental ex¬ 
pression for the nncleation landscape Vs{xc) that 
can easily be sampled in MC simulations and we 
show that the diffusion coefficients in the Fokker- 
Planck equation are a simple generalization of ex¬ 
pressions already availabl o®°i®^ for the fully trun¬ 
cated case = (ric)- Through our nonequi¬ 

librium statistical mechanics we obtain contri¬ 
butions to the nucleation barrier that have al¬ 
ready been fonnd throngh mesoscopic irreversible 
thermodynamics arguments concerning momen¬ 
tum corrections and by McCraw and Laviolette 
for energy fluctuations. In particular, our frame¬ 
work allows to rigorously define a nucleus tem¬ 
perature Tc that obeyes the non-Gaussian behav¬ 
ior derived by McCraw and Laviolette^ which 
has been confirmed in direct MD simulations of 
nncleatio n®®!®^ . We apply our approach to the 
most recent MD simulation predictions^ of Die¬ 
mand et al. and show that the lack of prediction 
of 3 — 5 orders of magnitude by the truncated the¬ 
ory can be eliminated. On the example of MFKT 
which uses the truncated state space we illustrate 
further that we can conveniently extend existing 
theories to richer state spaces thereby improving 
the quality of their predictions. 

In Sec. mi we develope the coarse-graining 
scheme for nncleation with extended state spaces 
from microscopic argnments. We show that the 
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nucleation landscape Vs is deeply related to the 
system entropy and that Vs is the fundamental 
static building block generating the nucleation 
dynamics on the coarse-grained scale. We discuss 
the nucleation landscape and the evolution equa¬ 
tion governing f{xc) in some detail and relate it 
to existing literature results. In Sec. [V] we apply 
our results to extend the state space of the MFKT 
of nucleation and discuss the nucleation of argon 
droplets in a metastable gas. Finally, in Sec. IVII 
we use our coarse-graining procedure to obtain 
the nucleation rate in metastable Lennard-Jones 
gas. We show that the standard coarse-graining 
procedure leads to erronoreous estimates as com¬ 
pared to brute-force MD simulations^ and that 
these can be corrected for by our theory. 


II. DETAILED DERIVATION OF COARSE-GRAINING 
PROCEDURE 

A. Relevant State Variables 

We aim to investigate nucleation in a metastable one 
component system. Prominent examples are nucleation 
of liquid drops in a gas phase or nucleation of gas bub¬ 
bles in a liquid phase. We characterize the state = 
(jiciPcUc) of the nuclei by means of their number of 
atoms ric, their total momentum Pc and their internal 
energy Uc- The incorporation of the total momentum 
allows for Brownian motion of the nuclei and the total 
energy as an additional variable is necessary to treat non- 
isothermal effects. The dynamics of a single nucleus is 
captured by the distribution function ft{xc) describing 
the probability that the nucleus has degrees of freedom 
Xc at time t\ the time-argument will subsequently be sup¬ 
pressed in the notation. The nuclei are assumed to evolve 
independently, so that the evolution of each nucleus does 
not affect other nuclei and the single-nucleus distribution 
function f{xc) is sufficient to capture the relevant macro¬ 
scale physics. We hence neglect the break-up of a nucleus 
into several nuclei or coalescence effects. We further char¬ 
acterize the full system in which nucleation takes place by 
its total energy Etot , total number of particles Ntot, total 
momentum Ptot and volume ftot so that its state is cap¬ 
tured by iCtot — (-^tot 1 .ftot, Ptot 5 fd;ot) • Clearly, as x^ot 
consists solely of conserved quantities, the state vector 
Xtot remains constant during the entire evolution. The 
state variables Xc (or respectively f{xc)) and a;tot are 
merely a splitting of degrees of freedom to capture the 
essential characteristics of the nucleus and its environ¬ 
ment. However, a;tot is not the state of the nucleus am¬ 
bient phase, but involves both, the nucleus and ambient 
phase degrees of freedom. In Sec. lIIEl we show that upon 
a suitable transformation, the variables a:tot and Xc can 
indeed be used to split the degrees of freedom into nu¬ 
cleus, Xc, and environment degrees of freedom, ccenv The 
splitting (a^tot, Xc) allows us to carry out the nonequilib¬ 


rium statistical mechanics more straightforward as a;tot 
is a conserved quantity. However, the state of the ambi¬ 
ent phase aJenv is not conserved as it exchanges particles, 
momentum and energy with the nucleus and as such, the 
nonequilibrium statistical mechanics is slightly more in¬ 
volved with the splitting (ajenv, Xc). 

The macroscopic state variables must be computable 
by exploiting specific microscopic expressions. On this 
more detailed level of description, the system is charac¬ 
terized by N atoms with mass m, position and mo¬ 
mentum Pi in a box of volume Vtot • All atoms build the 
microscopic phase space f = (ri,pi,..., rjv,PAr). In or¬ 
der to distinguish particles that are part of our nucleus of 
interest from the remaining ones we introduce the func¬ 
tion Hg^p) on r. We use the simple position dependent 
rule ng(p)(ri) = 1 if particle i is part of the nucleus and 0 
otherwise. Notice that in our notation ng(f)(ri) is seem¬ 
ingly only depending on r^, but this function of course 
depends on all position degrees of freedom of f. Popu¬ 
lar cluster criterions Hg^p^ of the kind are the Stillinger- 
cluster definition^, and its extension presented by ten 
Wolde et a h^^i^^ (see also Sec. I VI B]) . We obtain the sub¬ 
sequent microscopic definition of the macro-variables Xc 
and f{xc) characterizing the state of the nucleus: 

^Tic — n^j-p) (r^), 

i 

npe = Xlp*ne(f)(n), 

i 

™ jpi - 

““ ^ 2m 

i 

^f{xc) ~ ®c), 

where (j){ri) denotes the total interraction potential of 
particle i due to all other particles. The velocity 
of the nucleus is obtained from = Hp^. By 

averaging these quantities over phase space with the 
proper statistical ensemble, we recover the correspond¬ 
ing macrostate variables Xc- The variables = 

(HATtot, , Tlvtot) s-s the microscopic expressions 

for a^tot are introduced analogously. Using the shorthand 
notation H^,^ = (11™^ j Ilpe > n«c) '"^e also introduced the 
atomistic expression for the distribution function f{Xc) 
in Eq. ([5]). 




( 2 ) 


ne(r) 


B. Nonequilibrium Ensemble 

Under the asumption that the collection of variables 
(Ila^tot j 11/) captures all relevant physical processes on the 
macroscopic time-scale of interest, the nonequilibrium 
state of the system is characterized by the generalized 
canonical ensembl o'^^d^ : 

Pxtot,>‘c ~ ’ 7 ( 7^^(na;tot ~ 3;tot)x 

^ (^tot; j 


( 3 ) 
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where Z is the partition function of the mixed 
microcanonical-canonical ensemble Pxtot,Xa- particu¬ 
lar, the Lagrange multipliers Ac must be chosen to guar¬ 
antee for the proper macroscopic averages of the slow 
variables, i.e., f{xc) = (^f(xc)) Ac’ ’''^here the brack¬ 
ets denote an average over the entire phase space with 
respect to the ensemble Pxtot:>^c- We also have Xtot = 
(na:tot)pa.tot>^c ■ While the variables x^ot are treated mi- 
crocanonically since the total system is closed and the 
quantities aJtot rigorously conserved, the nucleus, as an 
open system, must be treated canonically through the 
use of Lagrange multipliers. 


C. Entropy 


We turn our attention to the entropy of the system 
described by {xtoti f{xc))- Subsequently, the entropy is 
shown to generate the irreversible dynamics of /. For the 
sake of clarity, we have summarized all different entropies 
involved in our analysis in Table ID 

We first derive a useful identity which relies on the 
observation that the nonequilibrium ensemble can be 
rewritten as 

Pxtat,^a — 7 / r^'^(na;tot ~ Xtot)e (4) 


and therefore the distribution / of the nucleus can be 
cast into the form 


f{Xc) {^fixc))p{xtot,Xc) 


Z(Xf^Ql , Xc) 
Z(xf^Qi^ , Ac) 


exp ( Ac(xc, Xtot)) ■ 


( 5 ) 


Here, Z is the microcanonical counterpart of Z, 


^(Xtot7 ^c) — / Xtot)f^(naj,, Xc), (6) 


and we associate it with the entropy S{xtot,Xc) = 
/cB log(.^(xtot, Xc)) of the system in state Xtot with the 
nucleus of interest having properties Xc. We obtain the 
entropy S'(xtot, /) of the syste m'^^d^ 


SjXtot,/) 

ks 


Iog(.^(xtot; Ac)) dXcAc(Xtot; ^tot)/*(®c) 


- / dXclog 


Ac(a3tot5®tot) 


f{Xc) 


- / dXc 


log (/(Xc)) - 


*5*(Xtot 7 


kB 


f{Xc) 


I ‘5'tot(^tot) 

(7) 


where S'tot(3Jtot) = log(^tot) is the entropy of the full 
metastable phase and contains its total partition function 


Ztot — dr Xtot) 



( 8 ) 


The contribution S'tot(a^tot) will turn out to be irrelevant 
in generating the final evolution equations. In the last 
line of Eq. ([5]) the integration is over the full phase space 
compatible with the system state Xtot. The expression 
for the entropy iS'(xtot, /) contains the entropic potential 


(ps{Xtot7 ^c) — (d(na;^ ®c))a:tot 

_ Ltot - Xc) (9) 



The entropic potential (jjs is related to the probability 
of occurence of a nucleus with properties Xc in a sys¬ 
tem described by Xtot- Upon studying the dynamics of 
the distribution / it will turn out that the entropic po¬ 
tential (jjs represents the nucleation landscape in which 
the diffusion process describing nucleation takes place, 
i.e., Us(xtot,®c) = -log(0s(xtot,a:c))- In the truncated 
case x^J^ = (ric) when the full system is isothermal at 
temperature T and Brownian motion of the nucleus is ir¬ 
relevant, the Gibbs-free energy of formation of a nucleus 
has been recognized2£“— to be /cbU times the logarithm 
of the probability that a nucleus with size ric appears 
in the metastable gas. Our expression Us(xtot,®c) = 
— log{(/>s(xtotj Xc)) hence represents a plausible extension 
to more nuclei degrees of freedom. In the general setup 
where the nucleus energy fluctuates, the Gibbs-free en¬ 
ergy of formation can no longer be regarded as the proper 
potential^^, but now. Us or respectively S'(xtot,iCc) in 
Eq. ([7]) as the entropy of formation takes this role. Eor a 
collection of the various forms of the nucleation landscape 
used throughout our work, see Table HIl 

The entropy of formation 5'(xtot,®c) is recognized to 
depend on the nucleus momentum Pc- To highlight this 
feature we show in App. IIX Al that in the limit of a suf¬ 
ficiently large nucleus ambient phase, 

3 

^(^totj^c) —‘8'(xtot; ric, ric) 2 log(nc) 

{^Vc{tIc 7 Pc) '^tot{Atot; Ptot)) 

— mUr - — -, 

c 27 - ’ 

( 10 ) 


where Utot = Ptot/{mN^ot) denotes the overall velocity 
of the system, — Pc/{mnc) the nucleus velocity and 
S'(xtot; n-c, Uc) represents the entropy of formation if only 
the variables {n^Uc) characterize the nucleus. 


D. Coarse Grained Evolution Equation 

Goarse-grained evolution equations can be split into a 
reversible part with the system energy as driving force 
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TABLE 1. Summary of relevant entropies. 
Expression Comments 


S{xtot, f{xc)) Total system entropy when state xtot and 
nucleus distribution f{xc) is known, see 
Eq. 0. Generates dynamics of f{xc) 
through Eq. (1571) . 

S'tot(a:tot) Entropy of full metastable phase, including 
the nucleus with arbitrary state. 

S{xtot,Xc) Entropy of system in state cctot containing 
nucleus in state Xc- 

S{xtot,nc,Uc.) Entropy of system in state xtot containing 
nucleus in state {nc,Uc)- Here, the momen¬ 
tum of the nucleus is unknown. See Eq. (nsj. 

S{xsnv,n.c,Uc) Entropy of system, where the system is split 
into a nucleus in state (ric, Uc) and its am¬ 
bient phase in state Xenv In this set of 
variables, the entropy can be decomposed: 
S{x 

envi^c^'l^c) — ‘S'env (:r^env) ~t“ Scij^c^ Uc). 

S'env(a^env) Entropy of nucleus ambient phase to dehne 
ambient temperature, pressure and chemical 
potential, see aSec. Ill El 


where the nucleation landscape Vs drives the dynamics 
and the elements of the diffusion tensor are 

D{Xtot,Xc) = — j , 

( 12 ) 

The average (•■•)(a;tot,a:c) configurations com¬ 

patible with the state Xtot of the total system and 
the initial state Xc of the nucleus. Here, AT-^j^Ha,^ = 
na;o('’'GK) — na,^( 0 ), where tqk is a time-scale on which 
the microscopic correlations drop off and the right hand- 
side in Eq. (HU reaches a plateau value. For the nu¬ 
cleation of liquid drops in a metastable gas, we rec¬ 
ognize TQK as the intermediate time-scale between the 
collision time of gas particles with the nucleus and 
the time-scale on which changes in the nucleus state 
are recognizable on the coarse-grained level of descrip¬ 
tion. We obtain therefore the tensorial generaliza¬ 
tion m of the expression for the diffusion coefficient 
(2rGK)”^ in Eq. (d]) already 

used in a variety of paper s ^i^ for the truncated state 
space x^J^ = (ric). The evolution equation (ED has been 
discussed in detail by Schweizer and Sagis^, but here 
we constructed microscopic expressions for its building 
blocks, D and Vs- 


Scij^c-, Uc) 


Entropy of nucleus with state x^ to d ehne E. Coarse Grained Evolution Equation in “Standard” 

nucleus temperature, see also Sec. (IHFI I. Variables 


Sc,K{nc,T) MFKT entropy of nucleus of size ric in am¬ 
bient phase at temperature T. 

Sc,K{nc,Uc) Extended MFKT entropy of nucleus in state 

{nc,Uc). 


and a part of irreversible origin involving the entropy 
as driving force^^i^. For our system the reversible time- 
evolution is absent since the only involved dynamics is the 
exchange of mass, momentum and energy of the nuclei 
with the ambient phase and these processes are known to 
be purely irreversible in homogeneous media. The funda¬ 
mental dynamics is hence encompassed in the irreversible 
contribution to the evolution equations. In App. lIXBl it 
is then shown that the evolution equations are recovered 
as: 


dxtot _ „ 
dt " ’ 


/ r dVs{xtot,Xc) \ 

\dxc dxc ) ’ 


( 11 ) 


In the literature nucleation is often described by de¬ 
composing the system into a nucleating phase and its 
ambient phas o^^’^^i^^ . We now aim to reformulate our 
theory in terms of the state variables (a^env,/), where 
the environment state Xem = (Aenv,Penv, Genv, Knv) is 
composed of the total number of particles iVenv) the in¬ 
ternal energy C/envj the momentum Penv and the volume 
I4nv of the nucleus ambient phase. Introducing the av¬ 
erage {A)f = f dXc f{xc)A of a coarse-grained variable 
A(xc, Xtot) with respect to the distribution function /, 
we can relate Xgnv to Xtot, i.e., 


I^env — I'tot ^ciXc)) f 

Pern — Ptot {Pc) f 

Uenv — Plot / ric “t” "Z 

\ 2mn, 
Aenv — Atot (^c) j- 


/ 


Pgnv 

‘luiN 0nv 


(13) 


The volume of the ambient phase I4nv is now a 
time-dependent variable since the volume of the nu¬ 
cleus Vc(xc) continuously changes during the nucle¬ 
ation process. In terms of the new variables the en¬ 
tropy S{Xf.oM,nc,Uc) with eliminated momentum degree 
of freedom, see also Eq. ED, can naturally be used 
to define the temperature T of the ambient phase as 
the derivative of S with respect to CAnv at fixed nu¬ 
cleus state, i.e., I/T = dS{xenv,'nc,Uc)/dUenv Anal¬ 
ogously, the chemical potential p, and the pressure 
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P are given by /r/T = —dS{xsnv,nc,Uc)dNenv and 
p/T = dS{x^^^,nc,Uc)/dVenv If we now assume that 
in the variables (a:env, Xc) the entropy splits into its 
environment and nucleus part, i.e., S'(a;env, nc, Wc) = 
*S'env(a;env) + SciucUc), where S'env(iCenv) Is a property 
of the ambient phase only and similarily the nucleus en¬ 
tropy ScijicUc) depends solely on the nucleus state, we 
find upon transforming the evolution equation dill) to the 
new variables. 


dt 


■ D{XenV) 

UXc 


Xc) 


f dVs{Xen^,Xa) \ 

\dXc dXc ) ' 

(14) 


Nucleation takes now place in the effective nucleation 
barrier 


TABLE II. Summary of relevant nucleation landscapes. 
Expression Comments 


Vs{xtot,Xc) Exact landscape valid for any system size, 
-log ((5(n,„, - 

Vs(a:env,®c) Exact landscape when the system entropy 
can be split into the nucleus entropy and am¬ 
bient phase entropy, see Eq. d. Used to 
extend MFKT to richer nucleus state space. 


Vs{xtot,Xc) Approximate Landscape valid when nucleus 
ambient phase is large, see Eq. (|27l). Used to 
study nucleation when nucleus environment 
is isothermal and isobaric. 


(^env; ^c) tic) 4“ rp rp pi 

{^cij^ci'Pc) ttenv (Nenv; T^env)) H , / \ 

+ nc --+ 2 ^b log(nc), 

(15) 

where ttenv = Penv/(mNenv) IS the velocity of the ambient 
phase. The expression dl5l) was already recovered in the 
nonequilibrium thermodynamics description of the nucle¬ 
ation process and has been discussed in full detail^. Our 
statistical mechanics approach however produces the ad¬ 
ditional logarithmic contribution in ric and the important 
observation is, that the introduction of the momentum Pc 
into the nucleus state space must be accompanied by this 
logarithmic term in the nucleation barrier. 


F. Entropic Potential in Thermodynamic Limit 


Sampling the entropic potential 4>s{xtotiXc) requires 
to calculate the probability of occurence of a nucleus with 
properties x^ in the microcanonical ensemble defined by 
the state x^^t- In the prominent case of a sufficiently large 
system ensemble averages become identical and instead 
of averaging over the microcanonical ensemble defined by 
a^tot we can sample over the isobaric-isothermal ensemble 
with fixed number of particles N, fixed pressure P and 
temperature T. By Galilean invariance we can set Ptot = 
0 without loss of generality. We define the corresponding 
system state by ®tot = {N, P,T, Ptot = 0). Averages of 
a phase-space variable Ua with respect to the isobaric- 
isothermal ensemble are given by 




/ dr 

J dr e~ ) 


(16) 


Vs(xtot,nc) Approximate truncated landscape valid 
when nucleus ambient phase is large. Used to 
study nucleation when nucleus environment 
is isothermal and isobaric. 

Vs(®tot, a^c®'^) Approximate landscape when nucleus ambi¬ 
ent phase is large and the nucleus described 
by its size and kinetic energy. See Eq. HSU. 
Leads to temperature fluctuation law of Mc- 
Graw and Laviolette. 


sample (/)s, such as the canonical or grandcanonical. Only 
when the system size is small, the rigorous exact version 
involving microcanonical averaging must be used. The 
momentum and position degrees of freedom can be sep¬ 
arated, so that the entropic potential can be sampled by 
MC simulations on the position of the particles. To this 
extent we integrate out the momentum degrees of free¬ 
dom and rewrite the entropic potential in a more suitable 
form: 


4^sixtot: Xc'j — dc^ {^5{\\-rLc rZc)d(np^ Pc) 


xS (n^pot -be), — Uc) 
^2 


x6 


n 


Pc 


2mn, 


POO 

oc / de'^F {nc,e'^,T)(j)s {xtot,nc,Uc-e'^), 
Jo 

(18) 


where 


For the entropic potential (jH]) we hence obtain the alter¬ 
native expression 

4‘s{XtotiXc) = ((5(na;^ — , (11) 

valid when the nucleus ambient phase is sufficiently large. 
By the same argument other ensembles can be used to 


n3/2r {Cv + 1) 



(19) 


and (ps (itot)ric,Uc — e(,) denotes the probability that in 
the isobaric-isothermal ensemble the nucleus of interest 
has size Uc and potential energy = Uc — e),. Here, 
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can be thought of an integration over all possible in¬ 
ternal kinetic energy contributions. We also introduced 
the atomistic expression Il^pot of the potential energy of 
the nucleus as the sum of the potential energies of all 
atoms of the nucleus, analogous to Eq. (ED, and the to¬ 
tal kinetic energy firkin, t of the nucleus as the sum of 
the kinetic energy of all its contained atoms. Further¬ 
more, Cv = fu-c ~ I is the heat-capacity of an ideal 
gas. The probability 4>s (aitot, rrc, rtc°*) can conveniently 
be sampled by means of MC simulations since it does not 
depend on the momentum degrees of freedom. 

To gain further insight, we approximate the probability 
4>s{,Xtot, ric, rtc°*) by the saddle point approximation ob¬ 
tained for each nucleus size ric in the potential minimum 
MP°*’°(nc) = min^^pot (/)5(aitot, ric, Me°*)- hr this minimum 
we can evaluate the second derivative of m-c, mP°*) 

with respect to which we denote by H{nc) and find: 


o^^JH{nc)4>s{xtot,nc) 


pot pot,0. 


X e 


— H{r] 


( 20 ) 


Here the factor of H{nc) is arising for normalization 
reasons. The entropic potential (/)s(itot, ric) of the trun¬ 
cated theory is recovered from 

= J dul°^ (j)s{xtot,nc,ul°^) 

= J dur^ - ne)5(n„po. - <°‘)) 

= (^(n«, - ric)) 

( 21 ) 


and with this we can define the nucleation landscape 
Vs(itot,Mc) = - log(())s(itot,Mc)) of the truncated the¬ 
ory. With the approximation (l20ll we can solve the inte¬ 
gral in Eq. (TTSl) analytically. A more elegant way is to 
rewrite 05 (®tot, ric) in the form 


05(^tot,^c) OC 2^2 


dei 


(Cy+ 1) 


V'i?(nc)0s(*t ot: 


-H 


T+Cv log 


( 22 ) 


We approximate the exponent in the integral to sec¬ 
ond order around the saddle point and we are then 
able to calculate the resulting Gaussian integral. The 
lower integration boundary can be shift to — oo since 
the main contribution of the integral comes from values 
> 0 so that the nucleation landscape Vs{xtot,Xc) = 


- log {<j)s{xtot,Xc)) is given by 


lf)?(^tot, ^c) — k5(^tot, Me) -t- 


Pc 


2mnckBT 


+ - log(nc) -I- log(r(Cy -I- 1)) + H- -1-^ 


-I- const (T) 


k^Tj 




+ o log 1 + 




y 




(23) 


that turns out to compare excellently with the true so¬ 
lution (1^ of the integral for nucleus sizes down to 
ric = 5, see also Fig O To gain further insight, we 





FIG. 2. Nucleation landscape Vs{xtot,Xc) of a liquid drop in 
a metastable argon gas for a nucleus with ric = 20,40 and 80 
particles and pc = 0 as a function of its internal energy Uc- 
The red line shows the approximate landscape from Eq. (I2ZD 
and the dotted blue line the exact landscape with integration 
of Eq. (I22II . Thermophysical data is taken from the sixth row 
in Table El 

can introduce a thermodynamic definition of the nu¬ 
cleus temperature Tc- We recall that Es(xtot,a:c) = 
— log(0s(:Etot(®tot), ®c)) and upon representing the nu¬ 
cleation landscape Vs in terms of the state of the ambient 
phase Xenv — (Aenv,Penv(— 0), Genv, f^nv), See ScC. Ill Ej 
we obtain 

dVs{Xtot,Xc) _ dVs{Xe,Xc) ^ dXe dVs{Xe,Xc) 

dxc dxc dxc dXf. , . 

_ dVs{Xe,Xc) _ dVs{Xe,Xc) 
dXc dXe 

Notice, that the variable set {x^, iCtot) is more convenient 
to evaluate the statistical mechanics, but the physical in¬ 
terpretation of Vs is more appealing in terms of (xc, Xg). 
When changing the nucleation landscape Vs{xe, Xc) with 
respect to the internal energy Uem of the ambient phase, 
while fixing the nucleus properties x^, this change must 
be equal to the temperature of the environment, such 
that dVs{xe,Xc)/dUe = —l/k^T, where the minus sign 
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is a result of our sign convention entering the definition 
of Vs through Vg = — log(0s), see also the form of Vg 
presented in Eq. (USD. Similarity, the nucleus temper¬ 
ature can be assigned as dVg{xe,Xc)/duc = —1/k^Tc. 
The change of variables law (IMl) implies 

dVg{xtot,Xc) _ _^ , 25 ) 

duc UbT UbTc 

and employing now the specific expression (I23p for the 
purpose of evaluating the left hand-side of Eq. (HSl) we 
find the caloric equation of state, 


Uc 


CvksTc + + 7^7 

kBH 




(26) 


and with this, the saddle point can be represented as 
e® = CyTc- To obtain Eq. (1^ we ignored the second 
last term in Eq. (051) since compared to the other terms 
it is of subleading order in ric- The internal energy of a 
nucleus is hence given by the kinetic contribution CyTc 
and a potential energy contribution involving the most 
probable potential energy V)? at a given nucleus size ric, 
and a correction for nonisothermal effects. When T = Tc, 
the nucleus is in the potential minimum of the landscape 
4>s{xtot,Xc) at fixed Uc and Pc- Plugging Eq. (051) and 
e® = CyTc back into the expression (051) for the entropic 
potential and expanding the logarithm of the T-function 
with Stirling’s approximation, we obtain the final form 
for the nucleation landscape 


P5(^tot7 ^c) -^c) T 


Pc 


2Tnc 


-a, 




+ 210 SK) 

-I- O f “ j + const(T). 


(27) 


Note, (j)g{x tot, ric) is the usual barrier to nucleation in 
the absence of Brownian-motion and non-isothermal ef¬ 
fects. In particular, 3/21og(nc) of the correction are due 
to the introduction of the momentum pc into the nucleus 
state space and l/21og(nc) through the Stirling approx¬ 
imation on the Gamma function r(C'y + 1). We here 
recover all corrections reported in the literature, includ¬ 
ing the momentum contribution obtained through en¬ 
tirely thermodynamic arguments by Reguera and Rubi^, 
and the temperature fluctuation law noticed by McGraw 
and Laviolette^. The additional temperature fluctuation 
contribution proportional to 1 /H has not been reported 
in the literature to our knowledge. In the subsequent 
studied model system of argon nucleation the contribu¬ 
tion proportional to 1/iJ effects the energetic landscape 
of Vg{xtot,Xc) dramatically and is hence relevant (see 
Eig. [2]). We postpone further discussion of the tempera¬ 
ture fluctuations to Sec. IlYl 


III. KINETIC ENERGY AS A STATE-VARIABLE 

In view of the various equivalent expressions for the 
nucleation landscape Vs, see Eq. (USD and Eq. ((27l) for in¬ 
stance, it is plausible to ask wether the kinetic energy of 
the nucleus can be used as a state variable instead of 
the nucleus momentum Pc, since the latter enters the en¬ 
tropic barrier only in the contribution p\l2mnc being the 
“systematic” kinetic energy of the nucleus. Suppose for 
the further considerations that Ptot = 0 by Galilean in¬ 
variance. If we use the modified set of nucleus state vari¬ 
ables x™°’^ = (nc, Uc) instead of Xc = (ric,Pc,Uc) 
we have to construct the entropic potential in terms of 
these new variables. The entropic potential (f>{xtot,Xc) 
is the probability that we encounter a nucleus with state 
Xc and therefore we find the clarrifying relationship 

(/)(a:tot,a:“°‘^) = J d^Pc S ^e^'"’® - ^ (l}{xtot,Xc). 

(28) 

This can be used to obtain the landscape 
Vs(a:tot,a:™°‘^) = - log((^(a;tot, which looks 

exactly like Vs{xtot,Xc), but with p\/2mnc replaced by 
the new independent state variable e^™’®. The modified 
set of state variables have the advantage that the 
state space is smaller upon catching the same relevant 
macro-physics. The intuition is that in a homogeneous 
medium it must be irrelevant in which direction the 
momentum Pc of the nucleus points and therefore the 
nucleation process can equally be described by \pc\, Pc 
or the kinetic energy e^™’®. 


IV. KINETIC ENERGY FLUCTUATIONS IN 
LITERATURE 

It has been convincingly demonstrated in MD 
simulations^!^ of nucleation that when only the prop¬ 
erties = {ric, e^™’ *) of a nucleus are measured, 

where now * is the total kinetic energy composed 
of internal and systematic contributions, the nucleus 
temperature Tc can be defined by Tc = 
and its fluctuations obey the law derived by McGraw 
and Laviolette, expIC^ [log (Tc/T) -I- (1 — Tc/T)]}, with¬ 
out the additional contribution involving H we obtained 
in Eq. (BID. We first remark, that the state space 
is a reduced version of Xc and therefore cannot capture 
entirely the rich physics contained in Xc- In Xc we have 
additional knowledge of the potential energy contribu¬ 
tions to the nucleus internal energy and it is then natural, 
that these contributions give rise to the more complicated 
caloric equation of state (l26l) rather than * = CyTc. 
Thus, in the reduced framework a;™'^ we measure ki¬ 
netic energy fluctuations only, and define the tempera¬ 
ture therein while in the richer theory Xc we measure the 
internal energy fluctuations composed of kinetic and po¬ 
tential energy fluctuations, and so it is not surprising that 
the fluctuation law contains additional contributions. 
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We apply a similar reduction of degrees of freedom as 
in Sec. m to show that in the reduced state space 
where only the size and total kinetic energy is measured 
in a simulation, the temperature indeed obeys the equa¬ 
tion of state ejj'”’ * = CyTc with the fluctuation law by 
McGraw and Laviolette. To this extent, notice that the 
integrand in Eq. m, denoted by (f){xtot,Xc, e^™’ *) rep¬ 
resents the probability that we observe a nucleus with 
property x^ and total kinetic energy Therefore 

= J d^Pcduc </>(itot, Xc, *). (29) 


Employing the approximation (uni the probability 
(j){xtot,Xc, *) has the form of the right hand side of 
Eq. (1^^ . but without the integration over e(,. Therefore 
the integrations in Eq. (j29p can be performed: 




^ kin, t / kin, t \ 

-l"C'vlog( j 

r{Cv + if 


so that the nucleation landscape Vs of the reduced theory 
satisfies up to irrelevant constants 


Vs(itot,3:c®‘^) =Vs(itot,nc) -l-log(r(Cv -I-1)) 


^kin, t 


fceT 


-I- Cv log 



(31) 


The temperature can be defined analogously to the pro¬ 
cedure in Sec. and we obtain e^'"’ * = CyTc and 
by employing the Stirling approximation finally the sug¬ 
gested fluctuation law by McGraw and Laviolette, 


Vs{xtot,x\ 


red 


) =Vs{xtot,nc) + - log(nc) 


-a 




(32) 


This shows that the kinetic energy of the nucleus is in¬ 
deed a valid measure to obtain the nucleus temperature 
if no information about its potential energy is available. 


V. APPLICATION: EXTENSION OF MFKT 

At the example of the MFKT of nucleationi^ we illus¬ 
trate that our work represents a powerful tool to extend 
existing nucleation theories to richer state spaces. 

In his worb^^, Kalikmanov derived a non-perturbative 
expression for the partition function Zc,K{nc,T) of a nu¬ 
cleus by microscopic considerations. In his treatment the 
ambient phase acts as a mean-field on the nucleus. His 
theory is valid for isothermal nucleation at temperature 
T and in particular he finds 

log(^c.K(nc,T)) = - ric 

-fcBK(ne,r)-l] 0 ,„iero(T) 

-f /CB l0g(VtotPsat('L)), 


where the chemical potential /Ltsat(iL) and particle den¬ 
sity PsatiT) at saturation are introduced. Here, n* can 
be interpreted as the number of particles located at the 
surface of the nucleus and 0micro is a “microscopic” sur¬ 
face tension. 

The entropy Sc,k{tIc,T) and internal energy udjicT) 
of the nucleus, see also Sec. Ill El are then given by 


Sc,¥.{nc, T) = kB log(Zc,K(nc, T)) -|- fcuT 

^2 dZc,K{nc,T) 


dZc^K{nc,T) 


dT 


Uc(nc,T) = kBT^ 


dT 


(34) 


by means of the usual definition of statistical mechanics. 


A. Enhancing State Space of MFKT 

We exploit a local equilibrium hypothesis to extend the 
entropy obtained by Kalikmanov to nonisothermal nucle¬ 
ation. To this extent, notice that the entropy S'c,K(nc, T) 
of a nucleus can be viewed as a property of the nucleus 
only and in this sense does not depend on the state of 
the environment. In isothermal nucleation the temper¬ 
ature of the nucleus Tc must coincide with the ambient 
temperature, so that T = Tc and we can equally write 
Sc,K = Sc,K{nciTc). Now the local equilibrium hypothe¬ 
sis enters and we postulate that a nucleus in nonisother¬ 
mal nucleation that can be assigned the temperature Tc 
has the same properties as a nucleus in isothermal nu¬ 
cleation that takes place at temperature T = Tc- Since 
the entropy, and analogously the internal energy of the 
nucleus are solely a property of the nucleus and do only 
depend on the nucleus state, we hence find the general¬ 
ization of the nucleus entropy to describe nonisothermal 
nucleation. 


Sc,K{nc, Tc) = kB log(Zc,K(nc, Tc)) + kBTc 


dZc,K{nc, Tc) 


dTc 


UciucTc) = kBT, 


2dZc,Kinc,Tc) 
dTc 


(35) 


where the nucleus temperature Tc is now related to its 
internal energy and mass by l/Tc{nc,Uc) = dS^^/duc- 
This then implies S'®))^(nc,M c) = S^'^^{nc,Tc{nc,Uc)). 


B. Analytic Solution for Nucleation Rate 

The entropy S^^(ric,Uc) does not depend on the mo¬ 
mentum variable Pc, so that we use the equations stated 
in Sec. Ill El to obtain the dynamics of extended MFKT. 
Our description then also takes naturally Brownian mo¬ 
tion into account. We take the diffusion tensor of the 
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particular forrn^^ 


D{nc) 


f(nc)l3x3 0 
V 0 ^ 


^ is not involved. This is however not true since the 
logarithmic normalization factor 3/2 log(n*) contained in 
AVs{xtot,x*) enters the prediction. 

For the prediction of the standard isothermal MFKT 
nucleation rate Jk we refer to the book of Kalikmanovi^ 
(Chapter 7.7 Steady State Nucleation Rate). 


and Z) is a 2 X 2 matrix with elements Dnc,uc 

and Du^,ua regulating the mass and heat exchange mech¬ 
anism between nucleus and ambient phase. The evolu¬ 
tion equations (I14[) of a nucleus in a metastable phase at 
temperature T and pressure P can be solved in the sad¬ 
dle point approximation^^i^i^ and give the prediction 
Jk, ext of our extended MFKT: 


Jk, ext —\ 


1 d^Vsixtot,Xc) 


Pg,i 


2tt dni 

1 I Amin I 
W y'IdetGI 


where 


G = 


e 


1 - 


^ 1 

c 

1 1 


V c c 

The coefficients are given by 


2 rp2 9uc{nc,T) 

O’,, — J 


dT 

1 /, d'^ns{T,nc) 

2 — L/microl,-t J 


( 7 : 


6.2 = 


drii 


(t:. 


C - ~CyT T 


,291og(psat(T)) 2 d'^ns{T,n) 


dT 


rpZ 


drindT 


(37) 


(38) 


-f 


dns{T, n) 


dur 


dOmicioiT) 

dT 


(39) 


and 


Amin = ^ (tr JZ - \/(tr^p^MdetH^ . (40) 


identified as the smallest eigenvalue of R = DG. All 
quantites in Eq. (1591) and Amin must be evaluated at 
the saddle point x* = = 0,m*) of the nucle¬ 

ation landscape Vs{xtot^Xc) as given by Eq. (ITKll . For 
the internal energy Uc{nc,Tc) this saddle point falls on 
u* = Uc{n*,T). Here, Al/s(®tot,iCc) denotes the height 
of the nucleation landscape and pg is the particle density 
of the metastable phase in which nuclei are formed. In 
the saddle point approximation (1571) the effects of Brow¬ 
nian motion are seemingly irrelevant since the coefficient 


C. Application to Nucleation of Argon 

As a specific model system we consider the nucle¬ 
ation of argon droplets in a supersaturated gas in the 
low-supersaturation regime studied by Diemand et al. 
through brute-force MD simulations^. This allows to 
compare our predictions directly with first hand exact re¬ 
sults on the nucleation rate. For the physical setups listed 
in Table HV] we have computed the isothermal MFKT nu¬ 
cleation rate Jk and the corrected rate Jk, ext and com¬ 
pared it to the MD simulations by Diemand. Since we 
compared our nucleation rate predictions to MD simu¬ 
lations, we did not take the real thermophysical data of 
argon, but data obtained through MD simulations^, see 
also Table IIIII 

In Fig. [5] we illustrated the predictions of the stan¬ 
dard isothermal MFKT and the extended version. The 
truncated theory seems to deviate systematically from 
the exact MD rates Jmd- Our theory succesfully elimi¬ 
nates this deviation indicating that we chose the proper 
set of variables to capture the nucleation process on the 
macroscale. Notice, that we did not include effects of an¬ 
gular momentum and shape fluctuations in our descrip¬ 
tion. As argued by Angelil et ali^ the effect of angular 
momentum seems only to be pronounced for small nuclei 
consisting of less than 16 atoms and all our critical sizes 
were beyond this size, so that angular momentum should 
not be used as an independent state variable for an au¬ 
tonomous description of the nucleation process in the low 
supersaturation regime. Angelil et al. recognized that the 
shape of large nuclei can deviate significantly from spher¬ 
ical objects. Using the shape of a nucleus as an indepen¬ 
dent fluctuating variable might offer additional insights 
into the relevance of this correction. 

Specifically, we evaluated the MFKT nucleation rate 
by resampling to the following relationships for which 
precise MD thermophysical data was available, see 
Table Ell We shortly summarize the parameters ap¬ 
pearing in the MFKT and closely follow the lines 
of its inventor-i^i2^. According to Kalikmanov, the 
microscopic surface tension is given by dmicro(J'c) = 
-\og{-B 2 {Tc)psa.t{Tc)/{kBTc)), where B 2 is the second 
virial coefficient of the gas. Introducing the macroscopic 
surface tension ■jsiTc), 'we can define the ratio uj{Tc) = 
ls{Tc) /{kBTc0nucro{Tc))- We further define C = 
and the coefficient A = Ni/Cj — ?>/A — 3/2, where Ni is 
the number of neighbours of particles in the nucleus. The 
number of neighbours in the nucleus is evaluated from 
the argon-liquid properties and can be approximated 
by Nx{Tc) = 5.5116?7(rc)2 + 6.138377(Tc) -b 1.275, where 
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ri(Tc) = %Pi{Tc)d\^. Here, dhs is given by the Barker- 
Henderson relation with the Weeks-Chandler-Anderson 
decomposition of the intermolecular potential, and rep¬ 
resents the effective hard sphere diameter in the the¬ 
ory of liquids, (ihs(Tc) = cliT bi/{a 2 Tc + a^)a with 
the Lennard-Jones parameters a = 3.822 x 10“^° m 
and e = 1.654 x 10“^^ J. The coefficients are given 
by ax = 0.56165-^, 02 = 0.60899-^, 03 = 0.92868-^ 
and bi = 0.9718-^. The mass of argon particles is 
m = 6.690 X 10“^® kg. The average number of par¬ 
ticles n*(nc,Tc) = a{nc,Tc)nc located at the surface 
of the nucleus is obtained by solving a cubic equation 
for the introduced coefficient a{nc,Tc). In particular, 
a = 3wC(l — + 3(DAC^(1 — For the 


J 

Logio — 

JMD 



FIG. 3. Comparison of argon nucleation rates J obtained 
from recent large-scale molecular dynamics simulations^, 
Jmd, with the original mean-field kinetic theory^^ (blue) and 
our extension (red). At each temperature T the rate compar¬ 
ison for various supersaturation ratios S is presented, see also 
Table HVl The solid lines are a guide to the eyes. 


kinetic coefficient Dn^,nc w® used the standard formula 
obtained from kinetic theory^i^S 


Dr, 


Psat 


^/2'Kink^T 


A(nc), 


(41) 


where A{nc) = 47 rr^ is the surface area of the nucleus 
and Tc = ( 3 pinc/ 47 r)^/^ its radius. The remaining co¬ 
efficients are given by predictions of Barret^i^, i.e., 
Drir^u^ = 2kBTDn^,n^ and Du^,u^ = Q{kBT)'^Dn^,n^ 
and are a result of kinetic considerations of the micro¬ 
scopic mass and energy exchange process between ambi¬ 
ent phase and nucleus. 


VI. APPLICATION: NUCLEATION IN 
LENNARD-JONES GAS 

We study the nucleation of liquid drops in a metastable 
Lennard-Jones gas in the low supersaturation regime 
considered by the latest large-scale MD simulations^ by 
Diemand et al.In order to compare our nucleation rate 
predictions obtained from our coarse-graining procedure 


TABLE III. Thermophysical data of argon^. 


k^T jt 

Psatcr^/e 

7so-^/e 



1.0 

2.55 X 10“^ 

0.453 

0.696 

-5.26 

0.8 

4.53 X 10“® 

0.863 

0.797 

-7.75 

0.6 

2.54 X lO-^ 

0.882 

0.882 

-12.9 

0.5 

2.54 X 10“® 

0.921 

0.921 

-18.15 


TABLE IV. Nucleation rate predictions of argon at various 
temperature T and supersaturations S. The initial density 
of the metastable gas is pg. Shown are the predictions Jmd 
of the brute-force MD simulations^ by Diemand et al. the 
predictions of the original MFKT— , Jk, and the predictions, 
JK,ext of our extension of MFKT. 

kBT/e S pgo^ logio(^) logio(^) logio(^) 


0.5 

72.8 

0.008 - 

-11.36 

± 

0.03 

-9.15 

-12.14 

0.5 

59.2 

0.007 - 

-13.01 

± 

0.02 

-10.41 

-13.44 

0.5 

48.7 

0.005 - 

-15.07 

± 

0.04 

-11.83 

-14.92 

0.6 

16.9 

0.008 - 

-11.75 

± 

0.01 

-9.53 

-12.42 

0.6 

15.6 

0.007 - 

-12.64 

± 

0.02 

-10.17 

-13.10 

0.6 

14.0 

0.006 - 

-14.39 

± 

0.05 

-11.16 

-14.15 

0.6 

12.0 

0.005 - 

-16.32 

± 

0.73 

-12.88 

-15.93 

0.8 

4.02 

0.030 

-9.08 : 

±1 

3.01 

-7.36 

-10.08 

0.8 

3.55 

0.025 - 

-11.70 

± 

0.01 

-8.87 

-11.70 

1.0 

1.66 

0.062 - 

-10.94 

± 

0.01 

-8.83 

-11.79 

1.0 

1.63 

0.060 - 

-11.90 

± 

0.07 

-9.38 

-12.36 

1.0 

1.60 

0.058 - 

-12.92 

± 

0.07 

-9.96 

-13.00 


with the MD simulations, we employ the same Lennard- 
Jones parameters as Diemand et al. and consider here 
nucleation in the limit of a large ambient phase at tem¬ 
perature T and pressure P. We aim to obtain the static 
building block (j)sixc,Xtot) and the dynamic information 
Dxc,xa from atomistic simulations. This then determines 
the coarse-grained model. We apply the following scheme 
to obtain predictions of the nucleation rate: 

1. Sample configurations C distributed accord¬ 
ing to the isothermal-isobaric ensemble in 
Eq. (USD through MC simulations to obtain 
4>s{xtot,'nc,u^°^), here approximated by Eq. (1201) . 
Construct the nucleation barrier Vs{xtot,Xc) in 
Eq. (HZD. 

2. Use the configurations C as input in short-time MD 
simulations to obtain Dx^^x^- 

3. Solve the Fokker-Planck equation m in the sta¬ 
tionary limit to obtain the nucleation rate predic- 
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tion. 

Subsequently we offer more details to each of the steps. 

A. Lennard-Jones Model System 

We simulate the interaction between atoms through 
the long-ranged Lennard-Jones potential, expressed at a 
separating distance r by 

ifO<r<re, 

0 otherwise. 

The cutoff radius is rc = 5 x tr. Typical values of tr and e 
for a monoatomic Lennard-Jones system, made of argon¬ 
like particles, are a = 3.822 x 10“^° m and e = 1.654 x 
10“^^ J, with a particle’s mass m = 6.690 x 10~^^ kg. 
We also introduce the typical time-scale t = on 

which the atomistic dynamics takes place. We performed 
both, the MC and MD simulations with iVtot = 1000 
particles in the isobaric-isothermal ensemble. We used 
periodic boundary conditions in all directions. With 1000 
particles system size effects are negligible and we benefit 
from the fact that the system is small enough so that 
only one nucleus is formed in the system. 

B. Clustering of Particles 

For our simulations we need to formulate the mi¬ 
croscopic cluster criterion ne(f) introduced in Sec. Ill Al 
to evaluate nucleus properties. We here rely on the 
extension^ of ten Wolde et al. of the most popular 
and simple scheme, the Stillinger-cluster definition where 
nearby atoms, within a distance smaller than a cutoff dis¬ 
tance Tc = 1.5(7 are considered part of the same cluster. 
Ten Wolde et al. further distinguish between particles 
with a liquid like environment and vapor like environ¬ 
ment. A particle has a liquid like environemnt if it has at 
least 5 neighbours as defined by the Stillinger-criterion. 
All liquid like particles which are connected are consid¬ 
ered part of the same nucleus. The scheme by ten Wolde 
et al. is computationally cheap and has been shown by 
Wedekind et ali^ to offer several advantages over the 
simple Stillinger definition. The developed nonequilib¬ 
rium statistical mechanics does not indicate which cluster 
criterion should be used, but once the cluster criterion is 
fixed, the formulation of our coarse-graining procedure is 
independent of it. 

C. MC Simulations 

To sample the entropic potential we exploit stan¬ 
dard MC simulations^ in the isobaric-isothermal ensem¬ 
ble with fixed number of particles A^, pressure P and 
temperature T. The formation of a nucleus in the low 


supersaturation regime of the Lennard-Jones gas is a 
rare-event and so we resample to an importance sampling 
technique to sample regions in phase space that contain 
a nucleus. We use the umbrella-sampling technique to 
bias the simulation with the bias potential 

n„,,bias = (n„^ - ricf (42) 

that favours the presence of a nucleus with size ric in 
the system. We use different umbrella windows corre¬ 
sponding to different choices of ric- In particular, we take 
umbrella windows from cluster size 10 to 140 usually in 
steps of 10. The bias strength is chosen as k = 0.25. The 
cluster size 140 is in all simulations beyond the critical 
size n*. In our MC simulations of the NPT ensemble, 
we attempt at each step a volume changing move with 
probability 1/A. If the volume changing move does not 
take place, each particle is moved with a probability of 
5/A in this step. The maximum distance in any direc¬ 
tion a particle can move is set to 0.4 (t. This ensures in all 
our simulations an average acceptance probability of vol¬ 
ume moves and particle moves of 50%. To sample each 
umbrella window we first initialize homogeneously dis¬ 
tributed particles in space so that their density coincides 
with the density of the metastable state at the temper¬ 
ature and supersaturation we are investigating, see also 
Table El Subsequently we perform an equilibration run 
of 10^ MC steps and then start averaging for 10® MC 
steps to obtain the entropic potential. The cluster crite- 

Fs 



FIG. 4. Nucleation landscape Vs(itot,nc) for argon at tem¬ 
perature T = 0.7e/fcB and supersaturation S = 6.11 near the 
critical size ric ~ 35. Dots are obtained from Monte-Carlo 
simulations and the different colors indicate different umbrella 
windows. The solid line is a fit with Eq. (H51) . 

rion of ten Wolde et al. allows for a nucleus larger than 
6 particles only. We however fit the data obtained for 
(j)s between Uc = 6 and 140 and assume this fit can be 
extrapolated down to Uc = 1. We found excellent fits of 









Systematic Coarse-Graining in Nucleation Theory 


14 


the nucleation landscape Vs{xtot,Xc) by Eq. (EZl) with 


Vs{xtot,nc) = - \og{S){nc - 1) + rji (ni - 


+ V2 


(ni - 


(43) 


where rji and 772 depend on the temperature T and super¬ 
saturation S. Furthermore, we fit the potential minimum 
MP°*’°(nc) = eoric + esuV^ and H{nc) = see also 

Table El Here, eo can be interpreted as the bulk energy 
per particle of a nucleus, and Cs as its corresponding sur¬ 
face energy. In Fig. 2] we illustrate a typical nucleation 
landscape Vs(£tot, n-c) near the critical cluster size recon¬ 
structed from the MC simulations. 


D. MD Simulations 


The atomic positional configurations harvested from 
MC simulations are used as initial conditions for our 
short-time molecular dynamics simulations in the micro- 
canonical ensembles. We can equally use the isobaric- 
isothermal ensemble at temperature T and pressure P 
to sample the short-time correlations of fluctuations de¬ 
termining the diffusion coefficient since on these short 
time-scales the macroscopic state of the system cannot 
be affected. Microcanonical simulations are more robust 
since there is no thermostat that can affect the nucleus 
state. If the isobaric-isothermal ensemble is used, the 
thermostats corresponding to T and P must be ensured 
to produce system relaxation times of the thermostated 
variables that drop off much slower than the time-scale on 
which relevant microscopic correlations drop off. Other¬ 
wise, the nucleus momentum and energy can be affected 
through the thermostat on the same time-scale as micro¬ 
scopic collisions of atoms of the ambient phase and the 
nucleus take place, resulting then in erroneous estimates 
of the diffusion tensor—. For simplicity we have sampled 
the initial momentum of all atoms in our MD simula¬ 
tions with a Gaussian distribution at the temperature 
T we are investigating. In this work we omit studying 
the energy dependence of the diffusion tensor and solely 
capture its size dependence. To sample the diffusion ten¬ 
sor at the critical cluster size for any given temperature 
and supersaturation regime, we take 32 configurations 
harvested from MC runs that contain a nucleus of the 
critical size and then run the MD simulation for a to¬ 
tal time of 25r around five times larger than the colli¬ 
sion time-scale tqk ~ 5r of gas particles with the nu¬ 
cleus. By averaging the dynamics of all 32 simulations 
we are able to obtain the required mean-square displace¬ 
ments m- The corresponding correlation functions 

are illustrated in Fig. [3 and 
within the collision time-scale a linear regime is recogniz¬ 
able from which the diffusion coefficient can be obtained. 
In Table El we have listed all obtained diffusion coefH- 


f^TQK Hj:,. ^tqk IIj:,.)(j:ioi.j:c) 
2«,2/3 



FIG. 5. Diffusion matrix normalized by for nucleation 
of argon at temperature T and supersaturation S = 6.11. 
The dots denote the mean-square displacement of ric with ric 
(red), between ric with Uc (blue), between Uc with Uc and 
with Pc, where Pc is the momentum in i-direction. The solid 
lines, the slope of which enters the prediction of the diffusion 
tensor, are fits to the linear regime of these correlations. 


cients with the convention 






3x3 


0 




0 1 

^ ^ricU 

\ 0 D''’° 0 "^'° / 

\ -^ricUc Uc,Uc/ 




(44) 


We also evaluated the diffusion coefficients for the tem¬ 
perature T = 0.7 and supersaturation ratio 6.11 for nu¬ 
cleus sizes distinct from the critical size and realized that 
they are well approximated by Eq. dm) for nuclei down 
to 15 atoms. The scaling (l4^ was also valid in the post- 
critical regime. We hence assume the validity of Eq. (l4^ 
for all nucleus sizes. In his wor k^^i^^ , Barret gives ex¬ 
plicit expressions for the coefficients and 

derived by kinetic theory arguments, see also Sec. IV GI 
We did not find his theoretical predictions to match with 
our computer experiments. But we found rather good 
agreement of his prediction of with ours. 

In evaluating the coefficients 

’^o^ice that the mi¬ 
croscopic process behind its contribution is the exchange 
of mass, momentum and energy of the nucleus with 
the ambient phase. We adopted the cluster criterion 
of ten Wolde et al. in our MG simulations, and within 
this criterion there are always particles located at the 
surface of the nucleus which are not part of the nucleus 
since they have less than 5 neighbours. These particles 
can however stay at the surface of the nucleus for a 
considerable time-scale much larger than the collision 
time-scale tqk and even rapidly oscillate between being 
part of the nucleus and part of the ambient phase. These 
oscillations lead to significant noise effects that make it 
difficult to evaluate the true diffusion coefficients, which 
by nature should only be affected by particles that truly 
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TABLE V. Fit of the entropic potential (1271) and diffusion tensor (1441) with conventions and parameters introduced in Sec. IVI Cl 
and Sec. IVI PI In eo and Cs as well as hi the error is smaller than the number of presented digits. The error in D^t has been 
discussed in Sec. IVIDl Errors in and are smaller than 0.05 and errors in smaller than 0.1. 


kuT/e 

S 


m 

V2 

eo/e 

es/e 

/iie^ 

D°t 

-‘-^nc,Uc^ 

^Uc,Uc^ 

C’°me 

0.8 

2.98 

0.0200 

7.81 ±0.04 

-9.73 ±0.25 

-4.8 

5.4 

6.8 

0.089 

0.72 

1.79 

0.77 

0.8 

3.11 

0.0211 

7.86 ± 0.04 

-8.99 ±0.12 

-4.7 

5.3 

6.6 

0.132 

0.71 

1.50 

0.62 

0.8 

3.33 

0.0230 

7.64 ± 0.05 

-7.73 ±0.24 

-4.7 

5.1 

6.6 

0.229 

0.68 

1.15 

0.50 

0.8 

3.55 

0.0250 

7.64 ±0.07 

-7.50 ±0.31 

-4.6 

4.9 

6.5 

0.314 

0.68 

0.88 

0.39 

0.7 

5.30 

0.0106 

10.62 ±0.07 

-12.19 ±0.23 

-4.9 

5.1 

7.5 

0.054 

0.75 

1.42 

0.80 

0.7 

6.11 

0.0125 

10.27 ±0.05 

-9.27 ±0.38 

-4.8 

4.9 

8.5 

0.145 

0.74 

1.16 

0.37 

0.6 

16.90 

0.0080 

12.21 ±0.06 

-7.63 ±0.70 

-5.3 

5.7 

4.8 

0.056 

0.54 

0.89 

0.49 

0.6 

15.60 

0.0073 

12.04 ± 0.07 

-7.26 ±0.21 

-5.3 

5.8 

5.1 

0.038 

0.61 

1.03 

0.63 

0.6 

14.00 

0.0065 

11.89 ±0.05 

-6.04 ±0.18 

-5.3 

5.8 

5.3 

0.032 

0.68 

1.20 

0.77 

0.6 

11.95 

0.005 

12.06 ± 0.08 

-7.19 ±0.12 

-5.3 

5.7 

6.4 

0.025 

0.74 

1.42 

0.81 

0.5 

39.37 0.0021 

16.34 ± 0.09 

-11.08 ±0.32 

-5.6 

6.3 

5.1 

0.013 

0.70 

1.00 

0.62 

0.5 

43.31 

0.0023 

16.57 ±0.10 

-11.94 ±0.28 

-5.6 

6.2 

5.2 

0.012 

0.67 

0.85 

0.56 

0.5 

48.70 

0.0026 

16.74 ± 0.09 

-12.53 ±0.23 

-5.6 

6.3 

5.4 

0.012 

0.58 

0.69 

0.48 


collide, are absorbed or emitted by the nucleus and do 
not fall in the undesirable class of oscillating surface 
particles, an inherent drawback of the ten Wolde cluster 
criterion. The contribution to the diffusion coefficient 
from the oscillating particles should dissappear if we 
average over sufficiently many simulations. The resulting 
noise effects are understood by the fact that the number 
of attempts with which the outer shell particles oscillate 
between being a nucleus property and a gas property 
outnumbers the true collision attempts from the ambient 
phase. We found that we can significantly improve the 
statistics of by using instead 

of the ten Wolde cluster criterion, the Stillinger-cluster 
criterion only so that now the oscillating particles in 
the ten Wolde cluster definition belong to the nucleus 
as well and thereby, the nucleus state Xc only changes 
through desirable collision processes. It is expected 
that this procedure results in an error in the obtained 
diffusion coefficients since by deciding that the outer 
layers of the nucleus is now also a nucleus property, 
the original nucleus state obtained from the ten Wolde 
definition is modified. This error only affects the 
prefactor in Eq. (Hill . The Stillinger criterion 

gives typically nucleus sizes nStiUmger ^^2 

larger than the ten-Wolde nucleus size j^tenWoide^ 

Eq. (liS we set = Dn^,n^l‘nJ'^. Suppose we ob¬ 
tained Dn^^ria in our simulation around the critical size 
J.J*, Stillinger 0^277,*. Wolde Stillinger Criterion, 

then and if we use the 

ten-Wolde definition, we must get the same coefficient 
Dn^^nc but D° = This means that 

the two values obtained for are distinct by a factor 
(^*.stiiiinger/^*.ten Woide)3/2 _ 2 , SO that the nucleation 


rate prediction that depends linearly on is affected 
by an error of a factor of two. In nucleation theory such 
errors are completely negligible. 


E. Rate-Estimation 

The rate prediction for the truncated theory using 
a;*''® = (uc) as state variable is given by the classical 
formula 


^trc — 



1 

27r 


d‘^Vs{xtot,nc) 

dnl 





(45) 


where n* is the critical size of the truncated nucle¬ 
ation landscape Vs{xtot,n*) and At4'(®tot,n*) the bar¬ 
rier height, see also Eq. (l43l) . In order to obtain the 
nucleation rate of the full theory Jext, we performed 
Brownian dynamics (BD) simulations of the full Fokker- 
Planck equation (HU and the Fokker-Planck equation of 
the truncated theory (HD to obtain the correction to the 
nucleation rate (H51) through nonisothermal effects and 
Brownian motion. For the BD simulation we closely fol¬ 
low earlier worb^. In Table lyll we tabulated the different 
obtained rates and corrections anon-iso and obm to Jtrc 
through nonisothermal and Brownian effects of the nu¬ 
cleus. We roughly estimated these correction by first per¬ 
forming the BD simulations in an isothermal setup, but 
with the correction for Brownian motion activated to ob¬ 
tain obm and subsequently we performed BD simulations 
incorporating all state variables and the additional cor¬ 
rection was attributed to anon-iso- We observed that in 
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TABLE VI. Nucleation rate predictions for a metastable Lennard-Jones gas at various temperatures T and supersaturations 
S. Jmd is the prediction of MD simulations^ (an asterix on top of the data indicates that the rate is obtained through the 
universal scaling relation^). Jtrc is the prediction of the standard technique using the truncated state space and Jext uses the 
full state space. Onon-iso and obm estimates the correction of the nucleation rate Jtrc due to nonisothermal effects and Brownian 
motion of the nucleus respectively. 


fcsT/e 

S 

logic (^) 

logic {^) 

logic {^) 

1 ^ — ( C^TinTi-iso ^ 

iogio ( aV i 

logic 

0.8 

2.98 

-16 ±0.5 

-11.43 ±0.29 

-16.08 ±0.61 

-1.0 

-3.7 

0.8 

3.11 

-15.53* 

-11.33 ±0.20 

-16.03 ±0.42 

-1.0 

-3.6 

0.8 

3.33 

-13.47 ±0.07 

-9.71 ±0.27 

-14.34 ±0.58 

-1.2 

-3.4 

0.8 

3.55 

-11.90 ±0.02 

-8.59 ±0.33 

-13.20 ±0.60 

-1.4 

-3.2 

0.7 

5.30 

-16.24* 

-11.89 ±0.31 

-16.40 ±0.59 

-1.2 

-3.3 

0.7 

6.11 

-13.28* 

-10.43 ±0.30 

-14.83 ±0.68 

-1.3 

-3.1 

0.6 

16.90 

-11.96 ±0.01 

-8.52 ±0.30 

-12.21 ±0.62 

-1.4 

-2.3 

0.6 

15.60 

-12.81 ±0.04 

-9.11 ±0.16 

-12.86 ±0.35 

-1.3 

-2.5 

0.6 

14.00 

-14.58 ±0.07 

-10.53 ±0.14 

-14.39 ±0.29 

-1.3 

-2.6 

0.6 

11.95 

-16.45 ±0.98 

-11.83 ±0.18 

-15.88 ±0.36 

-1.2 

-2.8 

0.5 

39.37 

-16.77* 

-11.49 ±0.22 

-15.43 ±0.48 

-1.5 

-2.4 

0.5 

43.31 

-15.89* 

-10.80 ±TaWe0.21 

-14.82 ±0.40 

-1.7 

-2.3 

0.5 

48.70 

-15.27 ±0.05 

-10.03 ±0.17 

-13.97 ±0.34 

-1.7 

-2.5 


J 

Logio — 
cr^ T 




Log (5) 

- 18 '. ,j 0.3 

1.90 1.95 2.00 2.05 2.10 2.15 2.20 2.25 (^-l) 

FIG. 6. Comparison of argon nucleation rate predictions of 
the standard Monte-Carlo technique (blue) with predictions 
of our theory (red) where the nuclei are described through an 
extended state space. Also shown are the results of brute- 
force molecular dynamics simulations^ (black). Dotted lines 
are a guide to the eyes. 


lowering the ambient temperature T the nonisothermal 
effects become more pronounced and comparable with 
the Brownian effects. In the high temperature regime 
however, Brownian effects dominate. The lower the tem¬ 
perature the more latent heat is exchanged between the 
nucleus and ambient phase so that we expect this behav¬ 
ior. In Fig [6] we compare the predictions Jext with the 
brute-force MD rates Jmd obtained by Diemand et alj^ 
and plotted it in a manner that highlights the intrinsic 
universal scaling behavior suggested by Tanaka et alj^, 


where the MD nucleation rate is recognized as a univer¬ 
sal function of \og{S)/{Tk/T — 1)^-^ with Tk = 1.3I2e/fcB 
the critical temperature. The rates obtained by the stan¬ 
dard truncated theory deviate by 3 — 5 orders of magni¬ 
tude from the exact MD rates whereas the rates obtained 
within our extended theory yields a quite excellent agree¬ 
ment. 


VII. CONCLUSION AND OUTLOOK 

We have derived the nonequilibrium statistical me¬ 
chanics of nucleation with an extended state space within 
the GENERIC coarse-graining technique. We have 
shown that the truncated state space where only the 
number of atoms in a nucleus is a relevant variable leads 
to systematically erronoreous prediction of the nucleation 
rate of 3 — 5 orders of magnitude for the prototypical 
studied nucleation of liquid drops in a supersaturated 
Lennard-Jones gas. When we extend the nuclei degree 
of freedom to involve in addition its internal energy and 
momentum, this lack in prediction can be significantly 
improved. For the example of the MFKT of nucleation 
we also showed that truncated theories can easily be ex¬ 
tended to a richer state space to improve predictions. 
Within our approach we were rigorously able to derive 
the correction of the nucleation landscape due to energy 
fluctuations and momentum fluctuations, and in both 
cases we were able to reproduce the up to date results. 
This includes the non-Gaussian temperature fluctuations 
inherent for small nuclei following the law of McGraw 
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and Laviolette, and the mesoscopic irreversible thermo¬ 
dynamics contribution of Reguera and Rubi for Brownian 
motion. The latter, as we found, must be corrected by 
a logarithmic contribution in the cluster size that arrises 
from normalization reasons, and can probabily not be ob¬ 
tained by bare thermodynamic argumentations as done 
by Reguera and Rubi. Within our statistical mechanics 
approach we were also able to state the proper definition 
of the nucleus temperature that leads to the temperature 
fluctuations predicted by McGraw and Laviolette. 

We predicted the nucleation rate in a metastable 
Lennard-Jones gas in the experimentally relevant low su¬ 
persaturation regime and found excellent agreement with 
recent large-scale brute-force MD simulations. As argued 
by Diemand et al. when the nuclei size is smaller than 
about 16 atoms, nuclei tend to spin so that the angu¬ 
lar momentum would enter as an additional important 
state variable. In the majority of our simulations the 
critical size was well beyond 16 so that these corrections 
are irrelevant. From this point of view it might be in¬ 
teresting to include additional state variables into the 
description such as the angular momentum or the shape 
of the nucleus. In this respect our coarse-graining proce¬ 
dure could also be extended to multicomponent system 
for more realistic applications. In our work we relied on 
a cluster definition to group the atoms as members of the 
gas phase or the nucleus of interest. Some work has been 
done avoiding the introduction of an artificial cluster cri¬ 
terion and we also expect useful extensions of our work 
in this directions, possibly following the lines of Kusaka 
et al^i^. 
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IX. APPENDIX 


where 


A. Truncated Entropy Splitting Eq. (IIOII 


Pxtot-Xc O'- ^(naiJtot 3^tot)<^(na;<; 3'c)i (50) 


Without restriction of generality, we can assume 
ftot = 0 emphasizing Galilean invariance. When the en¬ 
vironment is sufficiently large, the variables Xtot reach 
their thermodynamic limit and by the equivalence of 
ensembles we then have S{xtot,Xc) = log(Z(xtot, a^c)), 
where 

f "Etot 

Z{xtot,Xc)= / dT - Xc)e 

JVtat (46) 

— e ZPc — 0,yrc). 

The second line follows by straight forward calculation. 
This can be related to the truncated partition function 
Z{xtot, n-c, Uc) which finally defines the truncated entropy 
S(xtot,nc,Uc) by means of 

^ (^tot) ^c) — d Pc Z (a?tot: ^c) 

/'7n\3/2 

= Z{xtot,ric,Pa = 0,nc) (47) 

= (y) ^ Z{xtot,Xc). 

The required relationship (US than follows by exploiting 
Galilean invariance again to shift to the general case of 
arbitrary Ptot- 


is the fully microcanonical ensemble, the projection op¬ 
erator can be expressed alternatively as 

=(nA)(a!;tot,a5c=n,„c) 

S f 

+ (Tlictot ~ -^tot) • / '^■^c/(®c)(nA)(a;tot,a;c)- 

(51) 


The first part projects onto the nucleus variables 
while the second part onto the system state Xtot- Com¬ 
plementary, we define Q = 1 — P as the projection onto 
the fast variables. Denoting by iC the Liouville operator 
on f, we notice that iCWxtot = 0- The fact that aJtot 
are conserved quantites also implies on the macroscale 
dxtot/dt = 0. The projection operator technique then 
implie a^^d^ the time-evolution for /, 


df{xc) 

dt 


—{i£Ilf(^xc))p{ ■®tOt ) 


dXcMpxc),f{3i 


SS{xtot,f) 

Sf{x'c) 


(52) 


where 


2^ PTGK 

~T~ ^^(n/(a!;^)G(u)nj(3,^))p(3;j^j_A<,)- 

Jo 

(53) 


B. Detailed Derivation of Evolution Equation 


We use the GENERIC projection operator 
techniqu o^^d^ in order to find the fundamental time 
evolution of the distribution function / of the nucleus 
of interest. To this extend a projection operator V is 
introduced that separates the macroscopically relevant 
contribution of an arbitrary observable IIa on the phase 
space r from the irrelevant ones to be eliminated - this 
corresponds to the separation of slow, relevant from fast, 
irrelevant degrees of freedom when switching from the 
finer-level of description involving all N atoms to the 
coarse-grained level characterized by (a:tot, /)■ Following 
the standard definition, we define the projection onto 
the slow variables 

Pn^ =(nyl)p(a;tot,A,) 

+ (n •^tot 

dXc (nj(a,_,) — f{Xc)) ^ (n.4)p(cctot,Ae)- 

(48) 

Using 


(49) 


Here, n/(^„) = is the fast time-evolution of 

H/ and G{u) = is the time-evolution operator, 

the time evolution equation for /, Eq. (l5^ . splits into 
a reversible part and an irreversible part involving the 
Green-Kubo coefficients ) and driven by the 

total entropy S{xtot,f ) of the system. The time-scale 
tqk is an intermediate time-scale between microscopic 
collisions and the mesoscopic time-scale on which the nu¬ 
cleation process is described in a coarse-grained way. 

Using the auxiliary properties 

t£n /(.,) = - ■ i£Ux ^, (54) 

and (l50l) the reversible part of the time-evolution can be 
cast into 

= {i£Ilf(^xc)) pi •^'tot i^c) 

(55) 

We now turn our attention back to the projection op¬ 
erator P. This operator, or its counterpart Q act on 
i£Adx^ in the evolution equation. Therefore, we must un¬ 
derstand the object which also enters the 


dfixc) 

dt 


Pxtot,>^c^fiXc) ~ f {.Xc)Pxtot,Xc} 
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reversible part of the evolution equation. For homoge¬ 
neous nucleation this reversible part is known to vanish, 
so that we already expect 




(56) 


This implies ViCIVxc = 0 and QiCIixc = The 

argument for Eq. (1561) can be made more formal by rec¬ 
ognizing that the average over *£ 113 ;^ must be taken over 
all states compatible with the microcanonical ensemble 
Pxtot,xcJ that is over all states for which we encounter 
a nucleus with properties x^. But if a certain atom¬ 
istic configuration C = {ri, ...,rN,pi, ...ppf) contains a 
nucleus with properties Xc, then also the mirror image 
C = (—ri,—rjv,pi, ...pTv) obeys this requirement. It 
is important that the momentum variables must not be 
mirror imaged, since otherwise the momentum IIp^ of C 
and C do not match and therefore C would not fall into 
the class x^. In C all particles that would be about to 
join the nucleus would be about to leave it in C and 
vice verca. In the same manner one can argue that in 


C a momentum influx is transformed to a momentum 
outflux in C and vice verca. The same argumentation 
holds with the energy flux and consequently we have 
iCWxA^) = —iC^xA^') therefore the validity of 
Eq. diSl). 

In general, the irreversible part of the evolution equa¬ 
tion can be cast into the form^i^ 


df{xc) 

dt 


irr 




SS{XtotJ) 

’ SfixA 

6S{xtot,f) 


-JZ-f So, 


(57) 


where 

2 ^ rTGK 

M{xc) = I du{iCIlxi.G{u)iCIlxA p{xtot,xtot) 

Jo 

~ 7w Z (^TGKna!;c^TGKna;c)p(a!;tot,a!tot)- 
^KbTQK 

(58) 


We hence derived the fundamental nucleus evolution 
equation based on microscopic considerations. 







